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K-Optimal  Chaos  Ant  Colony  Algorithm 
and  Its  Application  on  Dynamic  Route 
Guidance  System 


Hai  Yang 


Abstract  Dynamic  route  guidance  system  is  an  important  part  of  the  intelligent 
transportation  system;  the  core  part  of  which  is  optimal  path  algorithm.  This  paper 
has  analyzed  the  main  influencing  factors  on  the  choice  of  optimal  path,  then 
provided  an  improved  K-optimal  chaos  ant  colony  algorithm  (K-CACA).  The 
road  impedance  factor  in  K-CACA  is  based  on  the  length,  crowdedness,  condition, 
and  traffic  load  of  the  road  sections.  The  optimizing  procedure  of  the  algorithm  is 
speeded  up  by  introducing  the  included  angle  threshold  of  direction.  The  chaos 
perturbation  effectively  refrains  the  algorithm  from  trapping  into  local  optima.  The 
results  of  simulation  experiment  show  that  K-CACA  is  effective  and  has  much 
higher  capacity  of  global  optimization  than  Dijkstra  algorithm  and  basic  ant  colony 
algorithm  for  optimal  route  choice. 


27.1  Introduction 

Dynamic  route  guidance  system  (DRGS)  is  an  important  part  of  the  intelligent 
transportation  system  [1].  It  offers  the  useful  optimal  route  guidance  information  to 
the  driver  according  to  starting  and  destination  point.  The  aim  of  DRGS  is  to 
achieve  the  goals  of  improving  traffic  system,  voiding  traffic  jam,  reducing  the 
vehicles’  travel  time,  and  realizing  the  rational  distribution  of  traffic  flow  on  each 
road  section  by  guiding  the  driver’s  travel  decision  [2].  So  the  core  content  of 
DRGS  is  the  detection  of  optimal  path  in  the  traffic  network. 

With  the  constant  enlargement  of  urban  road  network,  the  graph  theoretic 
algorithms  such  as  Dijkstra  algorithm  and  mathematical  programming  methods 
are  difficult  to  satisfy,  the  real-time  requirements  of  DRGS.  Initially  proposed  by 
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Marco  Dorigo  in  1991,  the  ant  colony  algorithm  (ACO)  was  aiming  to  search  for  an 
optimal  path  in  a  graph  based  on  the  behavior  of  ants  seeking  a  path  between  ant 
nest  and  a  source  of  food.  ACO  has  been  successfully  applied  to  solve  different 
combinatorial  optimization  problems  such  as  TSP  [3]  and  also  is  suitable  for  DRGS 
because  of  its  robustness,  positive  feedback,  and  distributed  computing. 


27.2  The  Influencing  Factors  of  the  Optimal  Path  Selection 
27.2.1  The  Length  of  Road  Section 

The  length  of  road  section  is  the  most  common  weight  during  the  procedure  of 
finding  the  optimal  path  because  the  length  is  easy  to  access  and  presents  the  travel 
cost  directly  [4].  But  the  optimal  path  does  not  only  mean  the  shortest  path  in 
dynamic  route  guidance  system.  The  travel  time,  crowdedness,  conditions,  and 
traffic  loadings  of  the  roads  also  affect  the  driver’s  decision  of  optimal  path. 


27.2.2  The  Travel  Time  and  Crowdedness 


The  travel  time  is  more  important  than  the  length  of  road  section  in  the  real  urban 
traffic  network  when  drivers  decide  which  path  they  should  go.  The  travel  time 
factor  is  dynamic  according  to  the  vehicle  number  per  unit  time.  So  we  can  compute 
the  crowdedness  of  road  by  the  travel  time  factor  [5]. 

The  crowdedness  of  road  refers  to  the  situation  that  the  vehicles  on  the  road 
cannot  run  at  the  normal  speed  so  that  both  the  travel  time  and  parking  delay 
become  longer.  The  crowdedness  factor  can  be  expressed  as  [T(ij)  —  T0(ij)]/ 
T0 (/,  y),  where  T(ij)  is  the  real  travel  time  of  road  section  (ij)  and  T0(ij)  is  the 
normal  travel  time  of  road  section  (ij)  without  jam.  In  general,  the  smaller  T(ij)  is, 
the  lower  the  crowdedness  factor  is,  and  vice  versa. 


27.2.3  Traffic  Load 

The  traffic  load  refers  to  the  number  of  vehicle  staying  on  some  road  section  for  a 
moment.  Particularly,  the  traffic  load  per  kilometer  presents  the  traffic  density  [6] . 
So  the  traffic  load  factor  presents  the  change  of  traffic  situation  over  time  and  its 
equation  of  state  is  as  follows: 


dx(ij)  (0 


dt 


HiJ)  (0  O(iJ)  (t) 


(27.1) 
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where  x^jjf)  is  the  vehicle  number  on  road  section  (ij)  at  time  t.  fij)(t)  is 
the  number  of  vehicles  coming  into  road  section  (ij)  per  unit  time.  0^^(t)  is  the 
number  of  vehicles  leaving  from  road  section  (ij)  per  unit  time. 


27.3  K-Optimal  Chaos  Ant  Colony  Algorithm 

The  urban  traffic  network  is  usually  abstracted  as  a  weighted,  directed  graph 
G  =  (KA,R),  with  node  set  V,  road  section  set  A,  and  R  =  {rtj\(ij)  ^  A};  rtj 
denotes  the  road  impedance  factor,  that  is,  heuristic  information  of  road 
sections  [7]. 


27.3.1  The  Selection  of  Subsequent  Node  Based 
on  the  Direction 

Because  the  dynamic  route  guidance  system  provides  optimal  path  information 
according  to  the  specific  starting  point  and  destination  point,  the  procedure  of 
finding  optimal  path  is  obviously  directional  [8].  When  the  ant  selects  next  node, 
the  probability  of  succeeding  in  finding  optimal  path  is  larger,  if  the  ant’s  moving 
direction  is  close  to  the  direction  of  destination  point.  Otherwise,  the  probability  of 
succeeding  in  finding  optimal  path  is  lower,  if  the  ant’s  moving  direction  is  far  away 
from  the  direction  of  destination  point. 

Let  vs  denote  the  starting  node  with  the  coordinate  (xs,  ys);  vd  denotes  the 
destination  node  with  the  coordinate  (xd,  yd);  Vi  denotes  the  current  node  where 
the  ant  stays  with  the  coordinate  (xb  yO,  and  Vj  denotes  some  node  adjacent  to  Vi 

with  the  coordinate  (xj,  yj).  Let  ViVd  denote  the  direction  from  current  node  Vi  to 
destination  node  vd,  called  the  approximate  destination  direction.  Thus,  the 

included  angle  6^  of  road  section  (vi9  Vj)  and  vivd  is  as  follows: 


Oij  =  arctan 


kd-  k 


next 


1  +  kd  -  k 


next 


(27.2) 


y  i  —  y.  _ » 

where  kd  = - -  denotes  the  slope  of  the  approximate  destination  direction  ViVd . 

xd  -  Xi 

y .  —  y. 

knext  —  — - -  denotes  the  slope  of  road  section  (\h  vd. 

Xj  —  Xi 

In  order  to  enlarge  the  probability  of  succeeding  in  finding  optimal  path,  we  can 
set  an  included  angle  threshold  cj)  .  The  subsequent  nodes  whose  included  angle  is 
greater  than  the  threshold  (j)  will  be  omitted.  The  convergence  speed  of  K-optimal 
chaos  ant  colony  algorithm  (K-CACA)  will  be  accelerated  by  decreasing  <^,  and  the 
solutions  diversity  will  be  enlarged  by  increasing 
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27.3.2  The  Heuristic  Information  of  Road  Section 


In  the  new  algorithm,  the  road  impedance  factor,  which  is  gained  by  using  weighted 
summation  of  the  length,  crowdedness,  condition,  and  traffic  load  of  the  road 
sections  [9],  presents  the  heuristic  information  of  road  sections. 

Let  rtj  denote  the  road  impedance  factor  of  road  section  (v^,  vj): 


COi  •  J~.  +  Cl>2  •  ^7  +  C£>3  ’  Qjj  +  •  Sij 


(27.3) 


where  Ly  denotes  the  length  of  road  section  (Vi,  Vj);  Cy  denotes  the  crowdedness  of 
(Vi,  Vj);  Qij  denotes  the  condition  of  (vi?  Vj);  Sy  denotes  the  traffic  load  of  (vis  Vj), 
and  &>k(k  =  1,  2,  3,  4)  is  weight  factor. 


27.3.3  Pheromone  Updating  Strategy 

The  pheromone  updating  strategy  of  the  new  algorithm  is  as  follow: 

L 

Tij(t+l)  =  (l-p)-Tij(t)+Y,^'ij  i=h-..,N  j=l,...,K  (27.4) 

1=1 


where  p  is  the  evaporation  rate  of  pheromone. 


27.3.4  Chaos  Selection  Strategy 

When  ant  colony  algorithm  initializes,  the  strength  of  pheromones  on  every  road 
section  is  equal  so  that  the  possibility  of  every  road  section  is  equal.  It  is  difficult  for 
ant  colony  to  find  an  optimization  path;  besides,  the  convergence  speed  of  algo¬ 
rithm  is  slow.  So  using  chaos  operator  is  necessary  because  it  can  increase  the 
searching  efficiency  by  the  random  and  ergodic  of  chaos  [10]. 

Take  the  logistics  mapping,  for  example,  iterative  formula  is  as  follows: 

Zi+I  =  p  ■  Zi  ■  (1  -  Zi),  i  =  0, 1,2,  •  •  •,  /<e(2,4]  (27.5) 

where  p  is  the  control  parameter,  with  the  domain  between  (2,  4].  When  p  =  4, 
0  <  <  1,  logistic  function  is  the  full  mapping  in  (0,  1),  which  is  at  a  totally  chaos 

status.  A  chaos  sequence  will  be  created  by  the  iteration  and  then  be  converted  to  a 
chaos  ergodic  parameter  when  solving  optimization  problems  in  space.  After 
introducing  chaos  perturbation,  the  new  status  transition  strategy  becomes 
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ri(0  =  < 


;,-(i  •  /,,)•" 


E  [4W]“-«-(1+zd 

0,  otherwise 


,  j  ^  allowedw  and  6tj  <  (j) 


j^allowedk 


(27.6) 


In  the  formula  above,  Og  is  the  included  angle  of  road  section  (vi9  Vj)  and  the 
approximate  destination  direction  ViVd .  is  the  included  angle  threshold. 


27.3.5  K-Optimal  Paths  Guidance 

The  K-optimal  chaos  ant  colony  algorithm  proposed  in  this  paper  provides  drivers 
with  K-optimal  paths,  which  are  almost  exactly  the  same.  Drivers  decide  which 
path  they  want  to  adopt  by  themselves,  according  to  the  real  conditions  of  roads  and 
the  travel  preference.  The  value  of  parameter  K  is  referred  to  the  scale  of  the  traffic 
network,  generally,  K  =  3.  The  newly  improved  algorithm  can  relax  the  conflict  of 
traffic  jam  due  to  that  many  drivers  choose  the  same  optimal  path. 


27.3.6  The  Realization  Steps  ofK-CACA 


The  realization  steps  of  K-CACA  are  proposed  as  follows: 

Step  1 :  Set  the  iteration  number  n  =  0,  initialize  the  pheromone  of  all  road  section, 
and  let  K  =  0. 

Step  2:  Put  the  starting  node  into  tabuk(s),  compute  the  road  impedance  factors  of 
every  road,  and  then  choose  one  subsequent  node  j  to  move  according  to  p\-  and  the 
included  angle  threshold,  finally  put  node  j  into  tabuk(s). 

Step  3:  Update  the  pheromone  of  roads  the  ants  passed  by.  Take  record  the  current 
optimal  path  and  let  K  =  K  +  1 . 

Step  4:  If  K  equals  to  the  specific  value,  then  stop  and  go  to  step  6. 

Step  5:  Adjust  the  initial  parameters  according  to  the  latest  traffic  situations,  run  this 
algorithm  again  in  order  to  obtain  new  optimal  paths.  Go  to  step  2. 

Step  6:  Output  K  optimal  paths. 
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27.4  Experimental  Results 

We  successfully  achieve  this  arithmetic  and  carry  out  simulation  experiments  using 
the  road  network  of  Jinan  City  shown  in  Fig.  27. 1  from  Google  Earth.  The  topology 
of  the  road  network  shown  in  Fig.  27.1  is  demonstrated  in  Fig.  27.2.  The  road 
impedance  factors  demonstrated  in  Table  27.1  are  obtained  by  calculating  the 
weighted  summation  of  the  length,  crowdedness,  condition,  and  traffic  load  of  the 
road  sections  at  8:30  in  the  morning.  Let  m  denote  the  number  of  ants  and  set 
m  =  20,  a  =  1,  P  =  3,  and  p  =  0.3,  and  use  logistic  mapping.  The  aim  of 


Fig.  27.1  The  road 
network  of  Jinan  City 
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Fig.  27.2  The  topology 
of  road  network 
of  Jinan  City 


Table  27.1  The  road  impedance  factor  of  road  sections 


Road  section 

rii 

Road  section 

rii 

Road  section 

rij 

Road  section 

rij 

(A,  B) 

4.96 

(B,  C) 

2.16 

(C,  D) 

2.02 

(D,  E) 

3.13 

(E,F) 

2.87 

(F.G) 

3.01 

(G,  H) 

3.08 

(H,  J) 

3.97 

(J,Z) 

3.99 

(Z,  W) 

4.01 

(A,  N) 

4.86 

(N,  O) 

3.22 

(0,  P) 

2.13 

(D,  O) 

3.99 

(P.E) 

3.91 

(P.O) 

2.71 

(P,  Q) 

2.21 

(E,  P) 

3.63 

(Q.F) 

3.70 

(Q.  R) 

1.86 

(N,  T) 

3.72 

(T,  Q) 

3.91 

(X,  W) 

3.77 

(R,  S) 

3.64 

(S,  V) 

3.71 

(S,  M) 

3.57 

(M,  L) 

2.28 

(H,  I) 

1.13 

(K,  L) 

1.97 

(F,Q) 

4.11 

(L,  W) 

3.11 

(V,  X) 

2.99 
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Table  27.2  The  results  of  three  algorithms 


Algorithm 

Optimal  path 

Distance  (km) 

Average  arrival  time 

K-CACA 

A-N-T-Q-R-S-M-L-W 

13.67 

22  min  and  7  s 

Dijkstra 

A-N-O-P-Q-R-S-M-L-W 

14.75 

26  min  and  3  s 

Basic  ACO 

A-B-C-D-E-F-G-H-J-Z-W 

14.14 

23  min  and  21s 

simulation  experiment  is  to  find  the  optimal  path  from  node  A  to  node  W  using 
K-CACA,  Dijkstra,  and  basic  ant  colony  algorithm  (ACO)  separately.  The  results 
are  shown  in  Table  27.2. 

The  results  of  simulation  experiment  show  that  K-CACA  is  effective  and  has 
much  higher  capacity  of  global  optimization  than  Dijkstra  algorithm  and  basic  ant 
colony  algorithm  for  optimal  path  choice. 


27.5  Conclusion 

In  this  article,  K-optimal  chaos  ant  colony  algorithm  is  proposed  for  finding  optimal 
path  in  dynamic  route  guidance  system.  The  road  impedance  factor  is  introduced 
into  K-CACA  based  on  the  length,  crowdedness,  condition,  and  traffic  load  of  the 
road  sections.  The  included  angle  threshold  both  accelerates  the  convergence  speed 
and  enlarges  the  solutions  diversity.  The  chaos  perturbation  effectively  refrains  the 
algorithm  from  trapping  into  local  optima.  The  experiment  results  show  K-CACA 
is  much  more  suitable  for  DRGS. 
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